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^Abstract. 

Q\ • Recently I proposed a new calculation scheme of a partition function of an immersion object 
1 using path integral method and theory of soliton (J.Phys.A (1998) 31 Mar). I applied the scheme 
£jto problem of elastica in two-dimensional space and Willmore surface in three dimensional space. 
|^!lfi this article, I will apply the scheme to elastica in three dimensional space as a more physical 
^jnodel in polymer science. Then orbit space of the nonlinear Schrodinger and complex modified 
Korteweg-de Vries equations can be regarded as the functional space of the partition function. 

> ■ 

i/-) . §1. Introduction 

O 

O . 

Elastica problem in two dimensional space M? has long history [1,2]. It is known that by 
^observing a shape of thin elastic beam, James Bernoulli named the shape elastica. It might be 
OYegarded as birth of the elastica problem and germination of the mathematical physics, including 
"^the elliptic function theory, mode analysis, nonlinear science, elliptic differential theory, algebraic 
• analysis and so on. The elastica in IR 2 [1,2] is defined as a curve with the Bernoulli-Euler 
^functional, 



E= ds fc 2 , (1-1) 



o 

> 

v-jwhere k is its curvature. 

•jjj . Recently I presented a new calculation scheme of a partition function of non-stretching elasticas 
M 2 under the condition preserving its local length [3] . The partition function is formally defined 



as 



DXe~^ dsk \ (1-2) 



where DX is the Feynman measure for an affine vector of a point of the elastica X and (3 is 
the inverse of temperature. Goldstein and Petrich discovered that the virtual motion of non- 
stretching curve obeys the modified Korteweg-de Vries (MKdV) equation, 

d t k + \k 2 8 s k + dlk = 0, (1-3) 

and its hierarchy [5,6]. Using the Goldstein- Petrich scheme, I found that the functional space of 
the partition function (1-2) are completely represented by the MKdV equation (1-3). In other 
words, the MKdV flows conserves the energy functional (1-1). The functional space (1-2) is 
classified by the solutions of the MKdV equation (1-3). 
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After that, I applied this method to the Willmore surface in three dimensional space M 3 
[4]. Instead of the MKdV equation, there appears the modified Novikov-Veselov equation which 
classifies the functional space of the partition function. 

In this article, I will investigate a partition function of an elastica in IR 3 with the energy 
functional 



where k is a complex curvature of the elastica in IR 3 . I will also require that the elastica does not 
stretch. 

Then the partition function of an elastica in IR 3 with the energy (1-4) can be also evaluated. 
Due to the non-stretching condition, instead of Goldstein-Petrich scheme of the MKdV hierarchy 
[3,5,6], the Langer-Perline scheme of the nonlinear Schrodinger (NLS) hierarchy and the complex 
MKdV (CMKdV) hierarchy appears in the calculation of the partition function [7,8]. 

Whereas the NLS equation is well known as the integrable equation and investigated well, 
the properties of the CMKdV equation is not sufficiently studied. According to the result of 
Mohammad and Can [10], the different version of the CMKdV equation does not pass the Painleve 
test [10]. In this article, I will also argue the properties of the CMKdV equation and the relation 
between the CMKdV and the NLS equations. 

On the other hand, the study of elastic chain model of a large polymer is current [11]. Accord- 
ing to recent review of a large polymer [11], statistical mechanics of a polymer model is closely 
connected with the mathematical science. Due to the complexity, investigation of its properties is 
not simple in general. However it sometimes can be exactly performed owing to deep symmetry 
[11]. In fact an exact partition function of elastic chain with the energy functional (1-4) was 
obtained by Saito et al. using the path integral [12]. However they paid no attention upon isom- 
etry condition as thermal fluctuation of the path integration even though they required isometry 
condition after all computations; they summed allover configuration space without isometry con- 
dition rather than over restricted functional space. It should be noted that the constraint does 
not commute with such evaluation of the partition function in general. 

Thus as another limit, it is of interest to investigate the partition function with the energy 
(1-4) under the isometry condition. One of purposes of this article is to investigate the partition 
function of a non-stretching space curve with the energy functional (1-4) as a polymer model. 

Furthermore, a space curve in IR 3 also interests us from the viewpoint of the string theory 
[15]. Grinevich and Schmidt investigated closed condition of a space curve obeying the NLS 
equation because a kind of its complexfication becomes a surface with Kahler metric [14] . Thus 
the problem is associated with the string theory [15]. (However as I mentioned in ref.[3], it should 
be noted that the elastica absolutely differs from a string in the string theory, even though it 
influences the theory [15].) Thus although it is not main purpose, another hidden purpose of this 
article is to investigate the moduli of non-stretching curve in M 3 by taking into the consideration 
of such relation as a generalization to the surface problem [4,14]. 

The organization of this article is as follows. In §2, I will evaluate the partition function of 
non-stretching elastica in IR 3 . Section 3 gives a discussion of the results. 

§2. Partition Function of Non-stretching Elastica in R 3 

I will denote by C a shape of an elastica (a real one-dimensional curve) immersed in three 




(1-4) 
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dimensional space R 3 and by X(s) = (X 1 , X 2 , X 3 ) its afiine vector, 



S 1 3 X(s) EC CR 3 , d s n X(s + L) = <9 s n X(s), (neN + {0}) (2-1) 

where L is the length of the elastica s is a parameter of the curve and N is natural number. I 
consider a closed polymer in R 3 ; its center axis is a space curve C. Here I will fix the metric of 
the curve C induced from the natural metric of R 3 ; 

ds = VdXdX. (2-2) 

As I stated in ref.[3], a reader should not confuse an elastica with a "string" in a string theory; 
they are absolutely different. 

There is the orthonormal system along C, (n , n 1? n 2 ) with fixing n as the tangent unit vector; 
n = <9 S X, where d s := d/ds. We make them, first, satisfy the Frenet-Serret relation [16], 

(2-3) 
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Here k is the curvature, r is the Frenet-Serret torsion and they are functions of only s. We rotate 
the orthonormal frame SO(2) fixing a := n so that we obtain (a , ai, a 2 ) [17-19], 














V a 2/ 



d s ^aij = | -«i | | ai | . (2-4) 
where k± := /ccos6>, Ki := ksin.9 and 

6(s) := f r(s')ds'. (2-5) 

J S 

For convenience, we introduce a complex curvature as 

k := m + ik 2 = ke w . (2-6) 
In this article, I will deal with a non-stretching elastica in R 3 with the energy functional 

E 



= [ ds \k\ 2 , (2-7) 
Jo 



which I will also call Bernoulli-Euler functional [3]. 

It is worth while noting that in general, there appear other potential terms in the energy func- 
tional for a general elastic rod. For example, there might appear elastic torsion term, stretching 
term and so on. An elastica is usually defined as a curve realized as a stationary point of an 
energy functional related to an elastic rod, at least, in the meaning of the classical mechanics. 
Hence the word "elastica" sometimes has ambiguity. Depending upon the potential term, its 
shape might belong to individual class. Thus reader should not confuse the word "elastica" with 
another one in another context. In this article, the word of "elastica" is meaning of a curve with 
the Bernoulli-Euler functional (2-7). 
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The elastica I deal with here is a model of a polymer which can freely rotate around its center 
axis but does not stretching and is forced by the potential (2-7). In other words, I assume that 
the force from the elastic torsion can be negligible but stretching can not. Furthermore, I will 
neglect the kinetic term of the elastica. Physically speaking, I will consider the polymers in the 
liquid whose temperature is determined and viscosity is very large. I also suppose that each 
polymer behaves independently and interaction among them are neglected. 

Let the elastica closed and preserve its local infinitesimal length for even thermal fluctuation; 
it does not stretch. Under the conditions, I will consider this partition function of the elastica 
given as [3], 

DXexp 1-/3 [ ds \k\ 2 } . (2-8) 



/ 
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Following the calculation scheme which I proposed in refs.[3,4], I will evaluate the partition 
function (2-8) under the non-stretching condition. 

However there is trivial affine symmetry of the centroid and direction of the elastica and the 
partition function naturally diverges [3]. For an affine transformation (translation and rotation 
g eS0(3)), X(s) — > Xo + (7X(s), (Xo and g are constants of s), the curvature k and the Bernoulli- 
Euler functional (2-7) does not change; this is a gauge freedom and the energy functional (2-7) 
has infinitely degenerate states. In the path integral method, I must sum over all possible states, 
Z includes the integration over M 3 and naturally diverges. As well as the arguments in refs.[3,4], 
I will regularize it, 

Zreg = Vol(Aff)' (2 " 9) 

where Vol(Aff) is the volume of the space related to the affine transformation. By this regular- 
ization, I can concentrate the classification of shapes of elastica. 

Next I will investigate the condition preserving local length even for the thermal fluctuation. I 
will expand the affine vector around the point which is an extremum point of the Bernoulli-Euler 
functional (2-7). I will call the point quasi-classical point according to the semi-classical method 
in path integral [3]. In the path integral, I must pay attention to the higher perturbations of e 
in order to obtain an exact result. Hence I will assume that X is parameterized by a parameter 
t. I will express a perturbed affine vector X around an extremum point X qc i in the partition 
function (2-9) as [3,4,7,8], 

X(s,t) :=e eat X qcl (M), ec>tX qcl = X qcl -X + 0(e 2 ). (2-10) 

with the relation 

<9 t X qc i = u a + itiai + u 2 a 2 , u a (L) = u a (0), (a = 0,1,2), (2-11) 

where -u's are real function of s and t. I will regard (2-11) as virtual dynamics of the curve 
describing the thermal fluctuation [3]. As in refs. [3,7,8], due to the isometry condition, I require 
[dt,d s ] = for X. Since ds qc i := a/ <9 s X qc i<9 s X qc ids, the isometry condition exactly preserves, 
ds = ds qc i. Here I will note that the deformation (2-10) generally contains non-trivial ones 
through u a (s) and the "equation of motion" (2-11). 

Let us compute the non-stretching condition [<9 t , <9 s ]X qc i = 0. I will introduce "velocities" 
(d t (f>i,d t (f>2) as 

d t a = <9*<9 s X qcl := (<9 t 2 )ai - (d t (pi)a 2 . (2-12) 
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and 

d s d t X qc i = (d s u - k±ui - K 2 -u 2 )a + {d 8 u\ + kiu )sli + (d s u 2 + K 2 -u )a 2 . (2-13) 

From the condition, I have the relation between dt4> c {4>c '■= 4>i + l 4>2) and a complex "velocity", 
u c := u\ + ru 2 , 

d t 4> c = i(« qc i«o + d s u c ) = i^dd' 1 Re((~K^u c ) + d s u c ) =: Q(u c ). (2-14) 

Here I use the notation K qc \ := k\ + ift 2 and I introduce the pseudo-differential operator d s in the 
meaning of 

d s U = Re(K^[U c ) = (K qc iU^+KaldU c )/2, 

u = d~ x Re{~K^Tiu c ) = J ds'Re(K qc i(s')u c (s')). (2-15) 

In order to find the connection between <p c and k, I will also investigate the fluctuation of a a 
(a = 1,2). By the virtual dynamics of ao, differentiation of a a (a = 1,2) by t must have the 
form, 

d t SLi = -<9t0 2 a o - v& 2 , d t & 2 = <9t</>ia + ua 2 , (2-16) 

where v means the rotation in the plane spanned by a a (a = 1, 2). By requirement of the isometry, 
the virtual dynamics of a a is constrained as [dt, d s ]a a = (a = 1, 2), 

-d s d t 8Li = {d s d t 4>2 - k 2 v )a + (#t02«i)ai + {d t (j>2^i + d s v)a 2 , 

-d s d t a 2 = -(d s d t (f>i - K 2 v)a - (d t (/) 1 K 1 )a 1 - {d t 4>iK 2 + d s v)a 2 , 
-d t d s a! = d t Kia + (Kid t <h)a-i ~ («i<9t02)a 2 , 

-d t d s a 2 = d t K 2 a + (« 2 <9 t 2 )ai - (« 2 <9 t 2 )a 2 . (2-17) 
Hence I have the relation [7,8], 

StKqcl = -Q(St^). (2-18) 

Accordingly I have the relation between <9t«; and complex velocity u c as the "equation of motion" 
of the deformation satisfied with the isometry condition [7,8], 

<9t« qc i = ~Q 2 (u c ). (2-19) 

I will remark that Q 2 is known as the recursion operator of the NLS and CMKdV equations. 
For this non-stretching deformation, the Bernoulli-Euler functional (2-7) changes as 

J \K\ 2 ds = J (|K q cl| 2 + e(7^d<9 t K qc l + Kq C ld t K^d) 

+ e 2 ((\d t K qc i\ 2 + K^c\d 2 K, qc i + Kqdd 2 !^) H )ds 

= J - e(/vT<3 2 (w ) + K qcl Q 2 (u )) 

+ e 2 ((\d t K qc i\ 2 +K~[d 2 K qc i + /c qc i^/qd) H )ds 

=: £ qc i + S^E qcl + S^E qcl + ■■■ . (2-20) 
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Since I wish to expand the complex curvature k around the extremum point in the functional 
space, I will require the extremum condition [3], 



S^E qcl = 0. (2-21) 

In this method, I will sum the weight function over all extremum points. Since they are 
extremum rather than stationary points, they need not be realized in zero temperature. 

Noting the relation 8 s uq = (7%Jw c + K qc iu^)/2 and above notices, supposed that Tt^Q 2 (u c ) + 
K qc \Q 2 (u c ) could be regarded as another function ~K^c[u' c + K qc \u' c of the variation of the normal 
direction in (2-15), I might find the relation 



/ 



/ ds(d s u ) = 0. (2-22) 



I supposed that the deformation is described by one parameter t. However there is no re- 
quirement that I must go along with only one parameter t to characterize this system. In the 
calculation of the partition function, one must sum up the weight function over events if the 
possibility of occurrence of the events can be considerable. I will search for all possible extremum 
points. 

Furthermore in a micro canonical system at energy E 0: the entropy S of the system is defined 
as S := log Z\e=e () and can be regarded as the logarithm of the volume of the functional space. 
From primitive consideration, the dimension of the functional space in the statistical physics is 
related to the degrees of freedom corresponding to E and the degrees of freedom of the elastica 
are not finite and its dimension need not one. 

Along the line of the arguments of ref. [3] , I will give up to express the thermal fluctuation using 
only one parameter t and I will introduce the sequence for mathematical times t := (ti,t3,ts, 
• • • , t2n+i, ' ' ' ) m this system so that (2-22) is satisfied. I will redefine the fluctuation (2-10) and 
introduce infinite parameters family, which can sometimes become finite set as I will show later, 

X Si = e (1 ^^^ s ' 2 »+ l9l 2»+iX qd 

= X qcl + 5t2n+id t2n+1 X qcl + O(l/0), (2-23) 

n=0 

where e was replaced with (l/v / ^)5t 2n+ i and d i2n+1 X qc \ is expressed as 

a t2 „ +1 X e = 4 n) a + ^ n) a 1+ 4 n) a 2 , = d^Re^u^), u^=Q 2n (u^). (2-24) 

The virtual equations of motion for the deformation are expressed as 

d i2n+1 n = Q 2n {u^). (2-25) 

Thus (2-25) represents the thermal fluctuation which conserves the local length. 

However it should be noted that there are two manifest symmetries in this system; one exhibits 
the symmetry of choice of the origin s and another is for the symmetry of U(l) phase of k; the 
later one is the same as the choice of the sq at the integration (2-5). For the transformation 
k(s) — > e lt /t(s — t), the partition function is invariant. 
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I require that the virtual motions must include such manifest symmetries 

%!«qcl = d s K qc i (2-26) 

and 

5ti«qcl = i«qci- (2-27) 

As in refs.[3,4], instead of the single deformation parameter, I will assign the infinite dimen- 
sional parameters in (2-23) to those which fulfill this requirement; t := (t, t) = (ti, £3, • • • , £1, £3, • • • ) 
In terms of these, I will investigate the moduli space of the partition function (2-9). In other 
words I will give a minimal set of the virtual equations of motion, which is satisfied with this 
physical requirement that the deformation contains the manifest symmetries (2-26) and (2-27), 

%„ +1 ^ci = (-Q 2 ) n (9 sKqcl ), d- t2n+1 K qcl = -Q 2 (d f2 „_ lKqcl ), (n = 1, 2, • • •)• (2-28) 

<9* 2 „ +1 «qci = (-Q 2 ) n Kci), d t2n+1 n qcl = -Q 2 (a t2 „„ lKqcl ), (n = 1,2,-..). (2-29) 

They are the CMKdV and the NLS hierarchies respectively. 

As stated in the introduction, the properties of the CMKdV equation is not well-known as far 
as I know. It has not been concluded that it is soliton equation yet. However even though it 
might not be integrable, properties of the CMKdV hierarchy and the CMKdV equation are very 
regular as I show as follows. 

As the NLS hierarchy, a solution of the n-th CMKdV equation, 

dt 2n+1 ^a - (-Q 2 r(d s K qcl ) = 0, 

is satisfied with the simultaneous equations by introducing unknown parameter fan-i, 

dt 2n+1 ^a - Q 2 (a f2 „_ lKqcl ) = 0, d hn _ lKqcl = (-Q 2 ) n-1 (<9 s K qc i). (2-31) 

This is a kind of Backlund transformation. Thus by ladder calculations, it can be proved that 
the solutions of the higher order equations belonging to the CMKdV hierarchy are also satisfied 
with the CMKdV equation, 

dfK+ -\K\ 2 d s K + d*K = 0. (2-32) 

In other words, the nontrivial deformation obeys the CMKdV equation as the soliton hierarchy 
does. (One might have a question why the ladder relation terminates at i = £3 rather than t = i±. 
From (2-26) i\ is determined as £~i = s + so and thus i\ is not an unknown parameter in the sense 
of (2-31). Thus (2-32) is a minimal non-trivial equation.) Further it is worth while noting that 
for intrinsically real initial condition 

k G e ia °!R for all s and constant cto for s and t, the CMKdV equation is reduced to the MKdV 
equation. Thus the solution space of the CMKdV equation has the completely integrable region 
in the meaning of the soliton theory. 

Similarly, I have the NLS equation as nontrivial deformation of the NLS hierarchy, 

\d t K+ -\k\ 2 k + 8 2 s k = 0. (2-33) 
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(2-30) 



For the NLS equation, this reduction can be naturally justified in the Jacobi variety of the 
hyperelliptic curve as a solution space [20]. 

Furthermore it is a very remarkable fact that for the variation of i obeying the CMKdV 
equation, the Bernoulli-Euler functional (2-7) is invariant, 

-d{ J ds \K(s,t,i)\ 2 = J dsd s (^\K\ 4 + (d 2 s R)K + Rd 2 s K - \d sK \ 2 ) =0, (2-34) 

as the NLS flows conserves the first integral, 

d t J ds |/<s,t,f)| 2 = -i J dsd s {{d s R)K-Rd s K) = 0. (2-35) 

Thus regardless of integrability of the CMKdV equation, it was clarified that its properties are 
very regular. 

Here I will comment upon the result of Mohammad and Can [10]. They investigated the 
"complex MKdV" equation and concluded that it is not a soliton equation. However their 
"complex MKdV equation" is expressed as 

dtK + \d s {\K\ 2 K) + din = 0, (2-36) 

which is a kind of " complexification" of the MKdV equation but differs from (2-32). Thus their 
result does not directly affect the studies on the integrability of our CMKdV equation (2-32). 

As I obtained the isometry deformation which includes the manifest symmetries (2-26) and 
(2-27), I will consider, here, the partition function (2-9). 

Since the CMKdV and NLS problems are initial value problems, for any regular shape of 
elastica satisfied with the boundary conditions, the "time" t and t developments of the curvature 
are uniquely determined . Furthermore noting that if one gives the real value k, k goes on real in 
the "time" t development of the CMKdV equation (2-32) whereas for the NLS equation (2-33) 
its "time" t development includes the complex value due to the pure imaginary in the first term 
in (2-33). Thus the "time" dt and dt are expected orthogonal in the moduli of the CMKdV and 
NLS equations. The "time" developments of both equations differ each other. In other words, 
for a given regular curve, there exist individual families of the solutions of the CMKdV (2-32) 
and NLS (2-33) equations which contain the given curve as an initial condition. Due to relations 
(2-34) and (2-35), during the motion of t and f, the Bernoulli-Euler functional (2-7) does not 
change its value. Hence the deformation parameter t and t draw the trajectories of the functional 
space which have the same value of the Bernoulli-Euler functional (2-7). 

In the case that I immersed an elastica in M 2 , the thermal fluctuation obeys the MKdV 
equation and there appears only one sort of hierarchy or the MKdV hierarchy [3] . In this article, 
the codimension of the immersion of the elastica in IR 3 is two while the former problem is one 
[3]. Accordingly it is natural that there appear twice degrees of freedom of the elastica in M 2 , t 
and i for the elastica in M 3 . 

Thus I can formally estimate the functional space for each functional value. By the "time" 
development of £3 and £3, I can classify the functional space of the partition function (2-9). In 
other words by investigating the moduli of the CMKdV and NLS equations which are satisfied 
with the boundary conditions, 



k(0) = «(L), X qc i(0) =X qcl (L), 
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(2-37) 



the measure of the functional integral d/i can be decomposed as, 



dfi = ^2dfiE- (2-38) 



E 



So I denote by the set of these trajectories which occupy the same energy E. 
Hence the partition function can be represented as 

Z Teg = I d/xexp(-/?£) = J2 e M~f3E) f dfi E = J2 e M-/3E)Vol(E E ), (2-39) 

E E 

where 

Vol(S s ) = I d/j E (2-40) 

is the volume of the trajectories S E . 

It is known that any solutions of the NLS equation (2-33) can be also expressed by the 
hyperelliptic function and its modulus agrees with the modulus of the hyperelliptic curves [13,20]. 
Grinevich and Schmidt studied the moduli of the NLS equation (2-33) whose corresponding space 
curve is satisfied with the boundary condition (2-37). For the NLS equation (2-33), there are 
infinite Jacobi varieties who have the same energy E in general. Thus it is expected that the 
CMKdV equation is connected among these Jacobi varieties induced from the NLS equation 
(2-33). 

As well as the arguments in ref.[3], even though I introduced the infinite dimensional co- 
ordinates t in (2-23), they are reduced to finite dimensional space, as the Jacobi variety of a 
hyperelliptic curve with finite dimension is embedded in the universal grassmannian manifold. 
Using the genus g of the hyperelliptic curves, I will evaluate the subspace and submesure of 
(E E ,d^ E ) as NLS part. The NLS part (S| LS ,d^g LS ) can be decomposed as (Sg LS ,d^ LS ) = 

For the case of a solution represented by the hyperelliptic function of genus g which is satisfied 
with (2-37), d(i E LS ^ is expressed as dt 3 A dt 5 A • ■ ■ A dt 2g -i. For each point, there is CMKdV 
flows. Even though it has not been confirmed that the trajectories of the CMKdV equation are 
linear and regarded as the vector space, it is clear that their cotangent space is flat and can be 
regarded as the vector space locally. Thus I can locally express the measure of d/x^ as 

d/4f } = dt 3 A dt 5 A • • • A dt 2g -i A di 3 A t 5 A • • • A di 2g -i- (2-41) 

Here I remove dti A dti in the measure because it exhibits trivial symmetries [3]. (2-41) is a 
subset of the infinite dimensional deformation parameters t in (2-23). Hence (2-40) becomes 

Vo\(E E ) = Vol(4 3) ) = E L d ^E- ( 2 - 42 ) 

9 9 ~v 

By exchanging the coordinate dt^ and dtj of multi-times t, the volume of is estimated by 
the unit of the elastica length L. Since the dimension of the Bernoulli-Euler functional E is the 
inverse of length and (3/ [length] is order unit, the multiple of the length can be interpreted as 
the multiple of the inverse temperature (3~ . Hence the sum of terms with different dimensional 
volume which appear in (2-39) can be regarded as expansion of power of (3. 
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§3. Discussion 



In this article, I gave a calculation scheme of the partition function of elasticas in R 3 in terms 
of solutions of the CMKdV equation (2-32) and NLS equation (2-33). Even though I could not 
give a concrete form of the partition function (2-9), I showed that its formal expansion is given 
by (2-39). As I thought that this scheme is based upon the soliton theory in refs.[3,4], I can 
not deny that it might be beyond the integrable system. In fact the CMKdV equation might be 
connected with the deformation of the Jacobi variety induced by the NLS equation (2-33). Hence 
I believe that this formulation might shed a new light upon the theories of the immersed object 
and its quantization (or evaluation of the partition function). 

Here I will mention the knot configuration. Since the NLS and CMKdV equations are initial 
value problems, the solution space includes any configurations of a space curve in M 3 . In other 
words, they also include any knot configurations and so I need not pay any attention upon the 
ambient isotopy [21]. In fact the trajectories of NLS equation classify space curves immersed 
in IR 3 rather than ones embedded in M 3 ; crossings are allowed and its topology disables us to 
distinguish such knot invariances or ambient isotopy. Since the knot configuration is physically 
discriminated by means of long range force such as the electromagnetic force and this theory in 
this article does not include such force, this notion can be physically interpreted. If one wishes to 
consider the knot configuration in this system, it might be related to the gauged NLS equation 
[22]. 

Next I will give two comments on the CMKdV equation. First, one might have a question 
why I need the CMKdV equation whereas the solution space of the NLS equation includes any 
configurations of a space curve in R 3 . I have been dealing with the measure of the functional 
space. An uncountable set of R becomes M 2 if the elements are measurable and one can define 
1R 2 topology in the set. In the similar meaning, I need CMKdV equation in order to introduce 
the natural measure in the functional space. 

The solutions of the NLS equation are described in terms of the hyperelliptic functions [13,20]. 
A hyperelliptic curve is embedded in a Jacobi variety. The trajectories of the NLS equation, 
(ti,t 3 , • • ■ , £23-1)) form the vector structure of the Jacobi variety. The NLS flow, which obeys 
the NLS equation (2-33), covers a subset of Jacobi variety. In each Jacobi variety, there exists 
compact subset as orbits of NLS flows. The individual Jacobi varieties are distinguished by 
points in the Siegel upper half space [20]. Since the CMKdV flows are perpendicular with the 
NLSE flows, the CMKdV flows might connect the different Jacobi varieties of solutions of the 
NLS equation. Thus I will conjecture, as the second comment upon the CMKdV equation, that 
the moduli of the CMKdV equation might be realized in the Siegel upper half space. It reminds 
me of the facts the theta function of the elliptic curve obeys the heat equation over the Siegel 
upper half plane, which is not integrable in the sense of the kinematic theory such as the soliton 
theory. The integrability of the soliton theory is associated with the time inversion symmetry 
and time translational symmetry, and the solutions are acted by a (continuous) group. On the 
other hand, the solution space of the heat equation is acted by only semi-group and thus it is not 
" integrable" in general. By complexfication of the heat equation, imaginary time heat equation, 
or Schrodinger equation is kinematic equation and integrable in the sense of the kinematic theory. 
Thus even though the CMKdV equation includes integrable solutions as kinematic region [23], 
I have a question regarding role of the CMKdV equation in Jacobi varieties of the hyperelliptic 
curves; is it in the framework of the integrable system? However in this stage, I cannot explicitly 
express the role of the CMKdV equation because there are few studies on the CMKdV equation. 
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I state that the properties of the CMKdV should be investigated. 

Finally I will comment upon the higher dimensional elastica problem, e.g., an elastica in n- 
dimensional space C C M. n . The codimension of the elastica becomes n — 1 and thus instead of 
t = (t, t), there appear (n — 1) sets of infinite dimensional parameters t = (t^\ t^ 2 \ • • • , t^ 71-1 )). 
As there appeared [/(l)-bundle in this article, they represent the (n — 2) -dimensional inner sphere 
of sphere bundle over the elastica C and the normal radius direction of C. Thus there is naturally a 
principal bundle over C. In other words, one can add the group structure over the equations. Thus 
the generalized MKdV equation naturally appears [8,24] and it is expected that my computation 
scheme of the partition function can be extended. 
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